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Abstract 

We study the correlations of two-dimensional dipolar excitons in coupled quantum wells with a 
dipole - dipole repulsive interaction. We show that at low concentrations, the Bose degeneracy 
of the excitons is accompanied by strong multi-particle correlations and the system behaves as a 
Bose liquid. At high concentration the particles interaction suppresses quantum coherence and the 
system behaves similar to a classical liquid down to a temperature lower than typical for a Bose 
gas. We evaluate the interaction energy per particle and the resulting blue shift of the exciton 
luminescence that is a direct tool to measure the correlations. This theory can apply to other 
systems of bosons with extended interaction. 
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The investigation of superfluidity of ^He and Bose condensation of alkali atoms stimulated 
the development of the theory of non-ideal Bose gas. The theory has been developed for 
contact and hard sphere interactions that correspond to the hard core repulsion between 
helium atoms and between alkali atoms. |j] During the last few decades a very 

intensive search of Bose condensation of excitons in semiconductor heterostructures resulted 
in the study of a new type of a Bose system, namely, a system of indirect excitons in coupled 
quantum wells in which electrons and holes are located in different wells, insert in FigIT] (see 

nri 

Refs.p, l6| and references therein). Such an exciton system is not only two-dimensional but 
also the interaction between excitons is essentially a dipole - dipole repulsion that extends 
to a much larger distance than the exciton radius or the separation between the wells. 

The purpose of the present letter is to demonstrate that the relatively large range of 
the dipole - dipole interaction compared to the contact interaction of bosons results in a 
dramatic difference of the system behavior. One of the main new features is that at any 
small concentration the quantum degeneracy of the gas with the reduction of its temperature 
is accompanied by a setting in of multi-particle correlations. As a result the degenerate 
system cannot be considered as a weakly non-ideal Bose gas but it is rather a Bose liquid. 
The other feature is that at higher concentration the repulsion between the bosons squeezes 
the wave function of each one of them and significantly reduces their overlap. As a result 
the temperature of setting of a quantum coherence is reduced compared to what is expected 
from a Bose gas. 

The importance of the long range repulsion has been recently demonstrated by Zimmer- 
mann and Schindler. [3] They calculated the blue shift of the exciton luminescence line which 
is essentially the average interaction energy per one particle. With the help of numerical 
calculations Zimmermann and Schindler found a strong pair correlation between excitons 
that was equivalent to a depletion region around each one of them. This depletion reduces 
the average interaction energy by about 10 times compared to usually used mean field value! 
This result is not only fundamental but also has a practical importance for experimentalists 
because the measurement of the blue shift is one of the main ways to extract the exciton 
concentration. 

In general, the reduction of the pair correlation function to zero at small distance has been 



noticed long ago 



8| and the crucial new point is the large length scale of this reduction 



0. 



In this letter we use analytical methods to study the correlation between excitons in the 
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whole temperature - concentration plane. In regions where the system can be considered 
as gas and only pair correlation is important we obtain analytical expressions for the pair 
correlation function and for the average interaction energy. At low temperatures and high 
concentrations where the system behaves as a liquid we use qualitative physical arguments 
to obtain reasonable estimates. 

The interaction energy between excitons in coupled quantum wells is 

where L is the separation between the centers of the wells and k is the dielectric constant. 
According to Ref. [lo| the attractive Van der Waals interaction is negligible in practically 
important values of L (L > 10 nm). If any correlation between excitons is neglected then the 
average number of excitons in area element dPr is ndPr where n is the exciton concentration 
and the average interaction energy is described by the "plate capacitor formula": 

i^.„. = j U{r)n£r = . (2) 

This result means that the main contribution to Eint comes from the interaction between 
excitons at distance of the order of L: [e^ / nL){mxL'^) = Tcne'^L/K. The same expression is 
obtained from quantum mechanical calculation in the mean field approximation. 

A necessary condition for the previous calculation as well as for the calculations that 
follow is nL"^ <^ 1. In the opposite case the electron - electron and hole - hole repulsion 
is stronger than the electron - hole attraction and it is hardly possible to expect a stable 
exciton phase. Practically, for excitons in typical coupled quantum well systems, this means 
that n < 10^^ cm"2. 

In reality, the correlation between excitons is not small. Two excitons with an energy of 
the relative motion E can come to each other only at the distance larger than ro{E) where 
[/(ro) = E. Therefore, Eq.([2]) is valid only if the region of strong correlation is small, i.e., 
ro -C L. The average exciton energy is of the order of temperature T and the last condition 
is equivalent to T ^ e^/nL. This condition can also be written as At ^ 7rV26L where 
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At = TiTi\j2/ MT is the thermal wavelength and b = Ti^k/Mc^ . The expression for b differs 
from the Bohr radius as = ti^n/me^ only by replacement of the electron - hole reduced mass 
m = memh/lme + rrih) by the exciton mass M = nif. + m^. In GaAs/AlGaAs structures 
the electron effective mass nif, = 0.067 (in units of the free electron mass) and the hole 



effective mass at the bottom of hhl level in a quantum well = 0.14. This gives M ^ 0.21 
and 6^3 nm which allows us to assume in further calculations that h <^ L. Therefore 



At ^ 7rv26L means also that \t <^ ttL which justifies classical description of the exciton - 
exciton interaction. 

While the validity of Eq.(l2]) requires T ^ j kL ^ 140 K at L ^ lOnm, the exciton 
binding energy e^, puts an upper limit to the exciton temperature. For L 10 nm e;, ~ 6 
meV = 70 K, i.e., if T ^ / kL more than 50% of excitons are ionized. There are also 
significant quantum corrections to U{r) due to overlap of the exciton wave functions because 
the exciton radius is not much smaller than 10 nm. 

These estimates demonstrate that the accuracy of Eq.(l2]) is quite limited. 

At T < I kL the minimal distance between excitons ro is larger than L and therefore 
it becomes the characteristic length scale of the exciton - exciton interaction. At this scale 
Eq-dID can be simplified and 



Uir)-^, ^o(T)^^^J . (3) 
At r ~ ro the exciton correlation is strong and 

Eint = n I U{r)g{r)d\ (4) 



where g{r) is the pair correlation function. 

We now introduce two important dimensionless quantities, the first is the 'interaction 
parameter': tq/At, which is a measure of the quantum character of the interaction, and 
the second is the gas parameter: nrg, which describes the "diluteness" of the exciton gas. 
When T becomes smaller than / kL the exciton wave length is still smaller than tq, i.e. 
tq/Xt ^ 1, and interaction between excitons can be considered with the help of classical 
mechanics. In this case g(r) = e~'^^'^^^'^ [see, e.g., Ref. 13[ Sec. 32] that results in 



Function JeIx) is shown in the insert in FigUl In the extreme case, T -C e'^/nL, it becomes: 

.2r2\ 2/3 

(6) 



/ 2 r2\ ^/■^ 

Eir^t = 27rr(4/3)n (-—j T^/^ = 27rr(4/3)nr2T . 



This result can be understood in the following way. Around each exciton there is a depletion 
region of radius ~ Tq. Without repulsion this region would contain nrg excitons with an 




T 



FIG. 1: Typical dependence Eint{T) for two different concentrations rii > n2. The upper insert 
sliows the exciton system. In the lower insert there is function fsix) that defines the interaction 
energy in Eq.([5]). 

average energy T. The energy necessary to force all of them out of the region is of the order 
of nr^T. 

Further reduction of temperature can violate one of the two conditions of validity of 
Eq.([S]). First, tq/Xt oc T^/^ can become smaller than unity, i.e., the exciton wavelength 
becomes larger than the potential length scale and the pair correlation has to be described 
with quantum mechanics. Second, gas parameter nrQ oc can become larger than 

unity. This makes multi-particle scattering and multi-particle correlations important, which 
is characteristic not for a gas but for a liquid. Which violation happens first depends on the 
concentration. 

If n ^ {h/2L'^Y then with decrease of temperature the exciton wavelength At becomes of 
the order of and then larger than tq while tq is still much smaller than the average distance 
between excitons (i.e., hTq -C 1). In this case the wave function ip{r) describing relative 
motion of two excitons penetrates under the repulsion barrier between them and the effective 
minimal distance is smaller than tq. Given the energy of relative motion the probability 
to find one exciton at the distance r from the other is e where {■ ■ ■)e means the 

average over all directions of the wave vector keeping E constant. For a non-degenerate 
exciton gas the probability density for an exciton to have energy E is (1/T)e~^/"^. That is 

9{r) = - mr)?)Ee-''l^dE (7) 

where 'il){r) has to be normalized in such a way that g{r)\r^oo = 1- At high temperatures, 
when Xt <^ tq it is possible to make use of the ip{r) semiclassical approximation. This leads 
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to g{r) = e~'^(^)/'^ and results in Eq.([5]) as expected. In the opposite case when Ay ^ tq (but 
still within the classical statistics limit, Xt ^ n^^/^) V^(r) falls exponentially at r < L'^/b 
and is a slow function of r at r > L'^/b. The interaction potential falls like U{r) oc and 
therefore the main contribution to the integral in Eq.(jl]) comes from the region r ~ L'^/b. 
The exact solution of the Schrodinger equation in this region gives 

^(0 = -T-TtIttit kJ^], (8) 



where k = V ME/h. The logarithmic dependence of this wave function on the energy is 
very weak and with logarithmic accuracy it is possible to put k ~ I/At- Then 



4n _2 / 2L 



and substitution of Eq. ([9]) in Eq. (jl]) results in 

^ Mln2(LV6AT) ■ ^ 
It is important to note that in all cases considered so far Eint is larger than the interaction 
energy at an average distance between excitons. This means that the main contribution to 
Eint comes from rear fluctuations of the interaction energy and these fluctuations are of 
the order or larger than the average value. It is intriguing that while at higher tempera- 
tures, where the interaction is classical, there is a clear temperature dependence of Eint, as 
can be seen from Eq.(l5]), at lower temperatures, quantum interactions yield an almost T 
independent Eint, quite similar to the mean field result, but with a different coefficient. 

The gas parameter in the above quantum case is different from classical nr"^ one. The 
exciton - exciton scattering crosssection is 

a= (11) 

2kln\kLyb) ^ ' 

(in 2D case it has units of length). In other words, a is of the order of the wavelength 
in spite the fact that the length scale of the potential is smaller than the wavelength 
(compare Ref . [l^ . Problem 7 to Sec. 132). Respectively the mean free path of excitons is 
I ~ \v? {E^ / b\T) / nXx ■ The gas condition in the quantum case corresponds to the absence 
of correlation between different scattering events which means that the wavelength has to 
be much smaller than the mean free path, i.e., n -C \y?{E^ /bXj) / ^'r- This inequality is 
identical with Eint/T <ti 1 and also with the condition of non- degeneracy with the accuracy 



of the logarithmic correction. Practically, the logarithm is not large compared to unity 
and hence the gas state of the exciton system is non-degenerate while the degeneracy is 
accompanied by strong interactions and multi-particle correlations between the excitons. This 
is one noticeable outcome of the above results. 

On the other hand, if n ^ (6/2L^)^ then a decrease of the temperature leads to a violation 
of the condition utq ^ 1 while ro ^ At- That is, quantum corrections are negligible and 
the system behaves as a classical liquid. The dimensional analysis gives 

E.nt = ^ n^/V ( ^ n^^A . (12) 



K \ kT 

where according to Eq.([6]) f{x) = 27rr(4/3)x~^/^ at x ^ 1. 

Now with a further growth of the concentration or a reduction of the temperature nr^ 
comes close to unity and a free motion of excitons between collisions becomes impossible 
because each of them is confined in between its neighbors. This classical picture is valid as 
long as the size of the confinement region is much larger than the exciton thermal wavelength. 
In a quantum language this situation corresponds to an excited state of the exciton in a well 
formed by its neighbors. The size of the well R ~ ra"^/^ and the motion is still semi-classical 
as long as <^ R- The energy at the bottom of the well ~ ze^L'^/R^ {z is the number 
of nearest neighbors) is of the same order as the depth of the well. Wells for different 
excitons are different, they are not static and sometimes some excitons overcome or tunnel 
across the surrounding barrier. But at nr"^ > 1 these rear occasions don't affect the order of 
magnitude of the estimates. In general, this picture is similar to a simple classical liquid and 
formation of the wells is the starting point of formation of a short range order characteristic 
for liquids|l5]. Further reduction of temperature brings particles to lower levels in the 
wells and makes the wells more stable. Stronger confinement of the wave functions of each 
exciton reduces their overlap. At nrg > 1 the energy of the bottom of the well is larger than 
temperature and it is the main part of the interaction energy. Eventually at nr^ ^ 1 most 
of the particles are at the ground state in their corresponding wells and rough estimates can 
be done assuming that a short range order has been formed and each exciton is inside a 
unit cell of a crystal. (We emphasize that we mean formation of short range order typical in 
liquids 1^ but not crystallization and formation of long range order |l^.) Then the potential 
energy of an exciton that is displaced by small distance r from the center of the unit cell is 



and R = [nsin(27r/2;)]~^/^ is the radius of the first coordinate circle. The values of the 
constants can be estimated assuming a crystal order within first few coordinate circles and 
random distribution outside. The dependence of the correlation radius on the order appears 

n 

quite weak and for hexagonal unit cell Ci ~ 10.5, C2 ~ 15 (compare Ref. |17l|). If T ~ 
{e^L?' / K)n^/'^ (i.e., ~ 1) then r R and there is no short range order in the system. 
However if T <^ (^^-^^Z K,)n^^'^ then the short range order does exist and r <^ R. The energy of 



the ground state in the potential (1131) above its bottom is ~ \ 2C2{h^/MR^){e^LyKR^) and 
quantization is important when temperature is of the order or smaller than this value, i.e., 
T < in^n/M) v/2C2(L%V2/^,). The order or magnitude of r is controlled by the temperature 
or quantization and in any case when the short range order does exist the size of a single 
exciton wave function is much smaller than the distance between excitons. 

The inequality r <^ R allows us to make two conclusions. First, the bottom of potential 
of Eq.( fT3i) gives a good estimate for the interaction energy Eint ~ Se^L^n^/^/zt. Second, it 
is possible to estimate the overlap of the wave functions of adjacent excitons. If L = 10 nm 
and n = 3 X 10^^ cm~^ the estimate according to the wave function in the harmonic potential 
of Eq. f|T3l) gives for the overlap ~ 0.2. The actual value is even smaller because Eq. f|T3|) is 
valid only for small values of r and the potential barrier at larger r is actually steeper. Due 
to the small wave function overlap the temperature at which the phase coherence or/and 
spin coherence 18|] in the exciton system is set in is reduced compared to its the expected 
value ~ fi^n/M. This is interesting as it points to a non-monotonic dependence of the 
quantum coherence onset temperature on the concentration, and is suggests that a lower 
density exciton system can become quantum coherent at higher temperatures than a higher 
density system, which is a-priori non intuitive. 

All estimates made above open the possibility to make up a general picture that demon- 
strates the role of particle correlations in a dipolar exciton system at the whole n — T plane. 
This picture is presented in Fig. [2l Correlation is not important and the mean field approx- 
imation is applicable only in region I. In region II the exciton system can be considered as a 
classical gas with strong pair correlations. The difference between region III and II is that in 
the former exciton - exciton scattering is described by quantum mechanics that changes the 
pair correlation function. Reduction of the temperature from region III to region IV leads to 
degeneracy of the exciton system. But simultaneously a strong many - particle correlation 
is set up. The system cannot be considered as gas, and the mean free path does not exist. 
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e/KL 



T = (e'L- 




FIG. 2: Different role of exciton - exciton correlation in the exciton system. In region I correlation 
is not important. II is the region of classical gas with strong pair correlation. In the region III the 
system is non-degenerate gas with quantum exciton - exciton scattering and strong pair correlation. 
IV is the region of quantum liquid. V is the region of classical liquid down to very low temperature. 

Rather surprising is the existence of region V where the system behaves as a classical liquid 
down to temperatures much lower than the degeneracy temperature in the case of contact 
interaction ~ h^n/M. The reason is that strong repulsion between excitons squeezes the 
wave function of each exciton to a size smaller than the average distance between excitons. 
In this region there is short range order and with a reduction of temperature the correlation 
radius grows, i.e., the system may crystallize. However, contrary to regular classical liquids, 
the attractive part of the exciton - exciton interaction is negligible lOj and it is likely that if 
a long range order is settled it is not as a result of a phase transition but as gradual growth 
of the correlation radius. 

It is necessary to note that the lines separating the different regions in Fig|2] do not 
correspond to any sharp transition. Crossing of one of the lines by changing temperature or 
concentration leads to a gradual change of the correlation between excitons. The study of 
the phase diagram is beyond the scope of this paper. 

The exciton interaction energy that produces the luminescence blue shift has different 
value and different dependence on T orn or both of them in different regions. This difference 
makes the blue shift a sensitive tool for measuring correlations in such an exciton system. 

We want to emphasize once again that the existence of region V is the result of the 
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significant tail of the exciton - exciton repulsion potential. In case of a short range potential 
(e.g., hard circles) an overlap of particle wave functions competes with the repulsion and the 
region of classical behavior does not exist. This latter case applies to excitons in one well 
where there is no dipolc - dipole repulsion and an increase of the concentration leads to a 
Mott transition but not to a classical liquid. 

In conclusion, 2D exciton gas with dipole - dipole repulsion never can be considered as 
degenerate Bose gas with weak interaction. At high concentration the repulsion suppresses 
all quantum effects. We evaluated the luminescence blue shift which is a sensitive tool for 
measuring of correlations in the system. 

B.L. appreciates discussions with M. Stern and A. L. Efros. 
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